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We discuss the role of electronic correlations in the iron-based superconductor LiFeAs by studying 
the effects on band structure, mass enhancements, and Fermi surface in the framework of density 
functional theory combined with dynamical mean field theory calculations. We conclude that LiFeAs 
shows characteristics of a moderately correlated metal and that the strength of correlations is mainly 
controlled by the value of the Hund's rule coupling J. The hole pockets of the Fermi surface show a 
distinctive change in form and size with implications for the nesting properties. Our calculations are 
in good agreement with recent angle- resolved photoemission spectroscopy and de Haas-van Alphen 
experiments. 
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I. INTRODUCTION 

High temperature iron-based superconductors have 
been intensively studied since their discovery four years 
agofD Among the various known iron pnictide classes, 
the 111 family comprising LiFeAs and LiFeP shows spe- 
cially interesting features compared to the other fami- 
lies. Whereas superconductivity in many iron pnictide 
compounds develops in the vicinity of a spin-density- 
wave (SDW) state upon doping or application of ex- 
ternal pressure, LiFeAs and LiFeP (and LaFePO from 
the 1111 family) are non-magnetic and superconductiv- 
ity evolves without additional doping or applied pres- 
sure. Of special relevance is LiFeAs where Tc 18 kP^ 
compared to 6 K"^ in LiFeP and ^ 4 in 

LaFePO. In the 1111 and 122 family compounds (with 
LaFeAsO and BaFe2As2 as typical examples), the SDW 
order is generally attributed to sizable nesting of the elec- 
tron and hole Fermi pockets.^ For LiFeAs the situation 
is not quite as clear: band structure calculations using 
density functional theory (DFT) predict an antiferromag- 
netic ground state with stripelike order as in the other 
pnictides, albeit in a shallow energy minimum compared 
to the non-magnetic state. ^ ^ In contrast, angle-resolved 
photoemission spectroscopy (ARPES) measurements re- 
port only poor nesting.^ In fact, recent neutr on sc attering 
measurements find strong SDW fiuctuation^^^l^ with an 
incommensurate vectoi^ slightly shifted from the com- 
mensurate order observed in the other iron pnictide su- 
perconductors. Also functional renormalization group 
calculation^^ predict SDW fluctuations to be the dom- 
inant instability. On the other side, recent de Haas- van 
Alphen (dHvA) experiments claim to be in good agree- 
ment with DFT regarding the topology of the Fermi sur- 
face.C^ 

It is of particular interest to identify the role of elec- 
tronic correlations in this context. Starting from band 
structure calculations within DFT, we include correla- 
tions at the level of the dynamical mean field theory 
(DMFT) and analyse their effect on the electronic struc- 
ture of LiFeAs. The band structure of LiFeAs features 



two shallow hole pockets around the Gamma point which 
generate a large density of states and it has been sug- 
gested that this is essential fo r the way superconductiv- 
ity emerges in this compound. These features of the 
electronic structure can also be expected to be rather sus- 
ceptible to changes induced by correlations. Thus, this 
paper aims to single out the effects of correlations on the 
Fermi surface and the low-energy properties of LiFeAs. 



II. METHODS AND INTERACTION 
PARAMETERS 

Our calculations were performed using an 
LDA+DMFT implementation which combines elec- 
tronic structure calculations in the full potential 
linearized augmented plane wave (FLAPW) framework 
with DMFT.^^ The electronic structure calculations 
were executed in WIEN2k^^, where the self-consistency 
cycle employed 1080 k points in the irreducible Brillouin 
zone, using the local density approximation^^ (LDA) 
for the exchange-correlation potential. We base our 
calculations on the experimental crystal structure as 
obtained from X-ray diffraction data^^ with space 
group PA/nmm. For comparison, we also performed 
calculations on the structure given in Ref. |2] for which 
we gave mass enhancements in Tab. [T[ 

For the construction of localized Wannier-like orbitals 
for DMFT, an energy window ranging from -5.5 eV 
to 2.85 eV was chosen, comprising the Fe and As 
4p bands. For the solution of the DMFT impurity 
problem, we employ paramagnetic calculations with the 
strong-coupling continuous-time quantum Monte Carlo 
methocP^as implemented in the ALPS cod^^^^^and con- 
sider only the density-density terms of the Hund's cou- 
pling; we used 1 x 10^ Monte Carlo sweeps throughout 
our calculations. 

For the interaction parameters, we use the definitions 
oiU = and J = (F^ + F^)/U in terms of Slater 
integrals^^ F^, and the fu lly localized limit (FLL) double 
counting correction.'^^'^ 
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There is considerable disagreement in the hterature 
about the size of the interaction parameters in the 
iron pnictides; in particular the Coulomb interaction U 
strongly depends on the estimation method, whereas J 
is only moderately reduced from its atomic value. Self- 
consistent GW determination yields rather large num- 
bers (e.g. U = 4.9 eV, J = 0.76 eV for BaFesAss"^, 
with lower values being reported by constrained LDA 
(e.g. = 3.1 eV, J = 0.91 eV for LaFeAsCPSi) and 
constrained random-phase approximation (cRPA) (e.g. 
U = 2.9 eV, J = 0.79 eV for LaFeAsO^^). For LiFeAs, 
interaction parameters obtained from cRPA have been 
reported in Ref. 26 for two low-energy models, one con- 
structed for the Fe 3d bands only, the other one for a 
manifold containing Fe 3d and As 4p states. The choice 
of the model affects the value of the interaction parame- 
ters in two ways: a model with more bands renders the 
associated Wannier functions more localized and thereby 
increases the matrix elements of the interaction. Also, 
since the interaction strength is derived as a partially 
screened Coulomb interaction where screening channels 
within the low-energy space are subtracted, the exclu- 
sion of more screening channels in a model with more 
bands increases the interactions. This is reflected by 
very different interaction parameters for the two mod- 
els, U = 2.45 eV and J = 0.61 eV for the d model, 
U = 4.95 eV and J = 0.87 eV for the dp model. 

However, as pointed out in Ref. [271 the appropriate 
model for our LDA+DMFT approach is a hybrid model 
where the Wannier functions are constructed from a dp 
model but only d-d transitions are excluded from the 
screening since we only treat the d states as correlated in 
our DMFT procedure. This means that the d model sys- 
tematically underestimates the interactions for our setup 
whereas the dp model systematically overestimates them. 
In light of these uncertainties we report in the following 
results for U = 4 eV, J = 0.8 eV and include a discus- 
sion about the sensitivity of our results to the choice of 
interaction parameters in Sec. |III[ 

In order to obtain real-frequency spectra from the 
imaginary time Monte Carlo data we performed analytic 
continuation of the self energy using the classic Maxi- 
mum Entropy method.^ To avoid uncertainties from the 
analytic continuation, the effective masses and Fermi sur- 
faces are directly infered from the self energy on the Mat- 
subara axis: the mass enhancements read 



^*/^LDA = 1 



(1) 



where the derivative is extracted by fitting a fourth-order 
polynomial to the data for the lowest six Matsubara fre- 
quencies.^^ The same polynomial is used for the deter- 
mination of the Fermi surfaces where we make use of 
the fact that the imaginary and real axis meet at zero, 
S(cj = 0) = ^{iuj = 0). 

Calculations were performed at a temperature T = 
72.5 K{f3 = 160 eV-i). 

In the following, orbital characters are labelled in a 



coordinate system which is 45° rotated with respect to 
the crystallographic axes, i.e. x and y point to nearest 
Fe neighbors in the Fe-As plane. 



III. RESULTS 

In Figs.[l]and[2l we compare the momentum-integrated 
and momentum-resolved spectral function for LiFeAs ob- 
tained within LDA+DMFT with their LDA counter- 
parts, namely the density of states (Fig. [T]) and the LDA 
band energies (Fig. [2]). Note that the LDA bands in Fig. [2] 
were renormalized by a factor of 2.17 corresponding to 
the orbit ally averaged value of the mass renormalization. 
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FIG. 1. (Color online) Orbital-resolved comparison between 
LDA density of states (red dotted lines) and the LDA+DMFT 
spectral function A{uj) (black full lines). The interaction pa- 
rameters used were [7 = 4 eV, J = 0.8 eV. 

The momentum-integrated spectral function A{uj) 
shows a bandwidth reduction but no substantial spectral 
weight transfer, i.e. no formation of Hubbard bands. The 
momentum-resolved spectral function A(k, cj) in Fig. [2] 
displays well-defined excitations around the Fermi level 
and stronger correlation-induced broadening of the states 
at higher binding energies. The broadening affects the 
states below the Fermi level more strongly where co- 
herent quasiparticles can be identified down to approx. 
0.3 eV below Ep. For the states above Ep^ the crossover 
to rather diffuse structures occurs at approx. 0.7 eV. 
On a quantitative level, at our temperature T = 72.5 K, 
the scattering rates (or, equivalently, inverse lifetimes) 
-ImS(zO+) are small, see Tab. |lj supporting the picture 
of well-defined, long-lived quasiparticles. The renormal- 
ized LDA bands give a good approximation only close to 
the Fermi level (the mass enhancement in Eq. ([T]) holds 
strictly only at = 0); states away from Ep are less 
renormalized. 
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FIG. 2. (Color online) Momentum-resolved spectral function 
A(k, uj) together with LDA bands. For comparison, the LDA 
band energies are divided by the orbitally averaged value of 
the mass renormalization. The interaction parameters are the 
same as in Fig. [l] 



For the given interaction parameters, the self energy 
and spectral function thus show the characteristics of a 
Fermi liquid state in a metal with moderate correlations, 
a picture which also has been promoted for the 1111 and 
122 fami ly of iron pnictides in a number of previous pub- 
lications .'^^'^^M^ Note that for multiorbital systems with 
sufficiently strong J, the absence of the rotationally in- 
variant Hund's coupling in the calculation {i.e. the con- 
sideration of the density-density terms in J only) can 
lead to qualitatively wrong results by suppressing coher- 
ence and driving the system from a Fermi liquid into a 
non-Fermi liquid state. ^^ ^^ This is not observed here, 
indicating that the restriction to density-density terms 
in the Hund's coupling is not detrimental. We also esti- 
mated the temperature below which coherent quasipar- 
ticles form by calculating = i-e. the param- 

agnetic local spin susceptibility at imaginary time point 
13/2. In a Fermi liquid = ~ T^^ and by studying 
the temperature dependence of = one finds 

that for U = 4 eV and J = 0.8 eV it takes on a constant 
value only at low temperatures below ~ 100 K. 

The mass enhancements as given in Tab. [T| exhibit 
pronounced orbital dependence, with stronger mass en- 
hancement in the t2g orbitals d^y and dxz/dyz- As can 
also be seen from Fig. TJ the bandwidth W of the t2g 
orbitals is smaller, leading to a larger ratio U/W and 
to increased correlations in these orbitals. Analysis of a 
low-energy iron d tight-binding model obtained by con- 
sidering the localized Wannier orbitals shows that the 
diagonal nearest neighbor hopping for the dxy orbital, 
tNN{xy^xy)^ almost vanishes as the direct hopping from 
the iron-iron overlap and the indirect hopping from the 
iron-pnictogen-iron overlap have opposite signs and al- 
most cancel. Also the diagonal hoppings to further iron 



neighbors for dxy are small; this contributes to the lo- 
calization of the dxy quasiparticles and a stronger mass 
enhancement than in the other orbitals. The table lists 



the mass enhancements for both investigated structures 
which show some quantitative yet not qualitative differ- 
ences. We checked that in particular Fermi surfaces are 
practically not affected, though; we therefore continue to 
give results only for the structure from Ref. il8| 



TABLE I. Orbital-resolved quasiparticle weights Z, mass en- 
hancements m*/mLDA, and scattering rates -ImS(iO+) for in- 
teraction parameters [7 = 4 eV, J = 0.8 eV. The first (sec- 
ond) number in each cell refers to calculations performed on 
the structure from Ref. 18 (Ref. 2). 



Orbital 



d 



'xz I yz 



Z 0.57/0.53 0.64/0.60 0.36/0.31 0.42/0.36 

mVmLDA 1.74/1.88 1.57/1.67 2.78/3.24 2.39/2.78 

-ImE(iO+) (meV) 0.1/0.3 -1.0/-0.7 2.4/5.2 1.7/3.8 

For a comparison with ARPES measurements. Fig. [3] 
shows some cuts of the momentum-resolved spectral 
function A{k.,uj) along the paths given in Fig. 2e in 
Ref. 9, Qualitatively, we find good agreement; quantita- 
tively, the mass enhancement extracted from these cuts 
in Ref. "9" is 3.1. This value should be compared to the 
mass enhancements of the orbitals that contribute most 
to the spectral weight at low energy. As can be seen 
from Fig. [l] these are the t2g orbitals (the Cg orbitals 
show a dip around the Fermi level) with calculated mass 
enhancements of 2.4-2.8 (2.8-3.2, respectively). Thus, 
mass enhancements are in good agreement, with ARPES 
pointing to moderately larger interactions. We will come 
back to this point further below. 




FIG. 3. (Color online) LDA+DMFT momentum-resolved 
spectral function A{k^uj) for LiFeAs along the paths in the 
Brillouin zone given in Ref.|9| Fig. 2e. Interaction parameters 
as in Fig. [l] 

In summary we consider LiFeAs a metal in an inter- 
mediate range of interactions without significant spectral 
weight transfer. Mass renormalizations are close to what 
has been calculated and measured in the 1111 and 122 
systems; compared to e.g. LaFeAsO^^, coherent quasi- 
particles seem to form a lower temperatures, though, 
with the spin susceptibility approaching Fermi liquid-like 
behavior only below ~ 100 K. 

We now turn our attention to the discussion of the 
effects of correlations on the Fermi surfaces of LiFeAs 
which have been experimentally accessed by ARPES and 
dHvA^ measurements. Figs. |4] and [s] show the Fermi 
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surfaces in the kz = and kz = tt plane obtained within 
LDA and LDA+DMFT. The pockets around {k^.ky) = 
(0,0) are hole pockets while the ones around ikx^ky) = 
(7r,7r) are electron pockets (compare Fig. [2|. The most 
prominent effects of correlations are the shrinking of the 
middle dxz/dyz hole pocket which takes on a butterfly 
shape at /c^ = 0, and the increase of the outer d^y pocket 
whereas the electron pockets almost don't change in size 
or form. This observation is in agreement with previous 
calculations^'' and would support ARPES results. This 
analysis shows that correlations tend to weaken -if not 
suppress- nesting in this material. 



the dHvA measurements, we show in Fig. |6] calculated 
dHvA frequencies with respect to magnetic field angle 
as reported in Fig. 2c of Ref. [131 The dHvA frequen- 
cies correspond to extremal pocket sizes (orbits) that are 
observed at a given angle 6 with respect to the kz axis. 
The labelling of the orbits follows Ref. 13: orbits 1, 2, 
and 3 refer to the inner, middle, and outer hole pocket, 
and orbits 4 and 5 to the outer and inner electron pocket 
(see Fig. |6] (a)). In order to define pocket sizes within 
LDA+DMFT (Fig. [6] (b)) in view of the finite broaden- 
ing induced by the correlations, we track the maximum 
of A{uj = 0) through the Brillouin zone. 
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FIG. 4. (Color online) Fermi surface for kz = 0. Color code 
for orbital characters: red: dxy, green: dxz, blue: dyz. Inter- 
action parameters as in Fig. [l] 
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FIG. 5. (Color online) Fermi surface for kz = tt. Color code 
and interaction parameters as in Fig. |4] 

For the discussion of the electron pockets, we de- 
scribe the Fermi surface in terms of an inner and outer 
pocket rather than by two crossed ellipses-like pockets of 
equal size. This is motivated by the fact that spin-orbit 
(SO) coupling lifts the degeneracy between the ellipses 
and splits the electron pockets into an inner and outer 
sheet .^^^^ Note, however, that no SO coupling is taken 
into account in the present calculation. As one can see 
from the comparison of Fig. |4] and Fig.jS) the thus defined 
outer pocket has strong kz dispersion whereas the inner 
sheet depends only weakly on kz. 

In order to facilitate a quantitative comparison with 



FIG. 6. (Color online) dHvA frequencies with respect 
to magnetic field angle obtained within (a) LDA and (b) 
LDA+DMFT. The orbits refer to extremal pocket sizes where 
the pockets are identified as discussed in the text. Interaction 
parameters as in Fig. [l] Triangles (pulsed field) and squares 
(dc field) are experimental data taken from Ref. 1 131 

Compared to the calculated dHvA frequencies in 
Ref. |13 small differences are already visible on the LDA 
level (Fig. [g] (a)), e.g. the minima of orbit 2 and orbit 
5 for small angles are not degenerate anymore. This is 
probably an effect of differences in the determination of 
the Fermi surface (e.g. due to effects of a finite k-mesh) 
and illustrates the high sensitivity of the orbits to to de- 
tails of the calculation. 

As already seen in Figs. [4] and [5j the effect of corre- 
lations on the Fermi surface manifests itself mainly in 
a shrinking of the middle hole pocket, and, in order to 
preserve the electron count, an increase of the outer hole 
pocket size. This is reflected merely by a shift downwards 
of orbit 2 and a shift upwards of orbit 3 in Fig. |6] (b); a 
change in the warping would indicate a change in the kz 
structure of the pocket which is not to be expected from 
the local, i.e. /c-independent interaction in (single-site) 
DMFT. 

Analyzing the curvature and the size of the orbits, the 
authors of Ref. 13 attribute the experimentally measured 
frequencies to the electron Fermi surface sheets, where 
the two higher frequencies are assigned to orbits 5b and 
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4a and the lowest frequency is suspected to originate 
from orbit 5a!^ Our results support this interpretation: 
whereas the orbits 2a/5a and 2b/4a are of similar size in 
the LDA calculation, the correlations affect mainly the 
hole pockets and lift this (near-) degeneracy. As a re- 
sult, the electron orbits 2a and 2b are unlikely to give 
rise to the measured frequencies as their sizes are rather 
different from the measured data. This offers a reconcili- 
ation of the dHvA and ARPES experiments: the shrunk 
middle hole pocket is only seen in ARPES which finds 
a correlated metal with poor nesting together with siz- 
able mass renormalization. In contrast, the dHvA mea- 
surements resolve the (lighter) electron pocket sizes in 
LiFeAs that almost do not change under inclusion of cor- 
relation and therefore report good agreement with LDA. 
The large mass renormalizations (up to ~5) that are also 
measured in Ref. 13 suggest -even under consideration of 
a non-negligible electron-phonon contribution- a scenario 
of important electronic correlations which are correctly 
accounted for within LDA+DMFT. 

In LiFeP, which was also investigated in the same work 
Ref. [131 the situation is different. Whereas the middle 
hole pocket in LiFeAs is particulary shallow with a band 
energy of 64 meV above Ep dit the F point, the pocket in 
LiFeP is roughly of the same size as in LiFeAs, however 
the corresponding band energy at F is 155 meV. Hence, 
the dispersion in LiFeP is considerably stronger, render- 
ing the pocket less susceptible to band shifts as induced 
by the real part of the selfenergy, i.e. less susceptible to 
correlations. Consequently, the dHvA frequencies remain 
rather unchanged upon inclusion of correlations. 



Sensitivity to interaction parameters 

We analyze the sensitivity of our results to our choice 
of interaction parameters by applying some variation to 
U and J while keeping the respective other parameter 
fixed. 
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FIG. 7. (Color online) Sensitivity of effective masses 
m*/mLDA with respect to changes in the interaction param- 
eters. 

Fig. [T] shows the evolution of the mass enhancements 
m*/mLDA with the interaction parameters. A moderate 
dependence on U and a very strong dependence on J are 



observed (note that the applied variations of U and J are 
different in size). Whereas a change in U affects all or- 
bitals roughly equally, a change in J leads to an immense 
mass enhancement particulary of the t2g orbitals. 

The decisive role of the Hund's coupling for the phys- 
ical properties of the iron pnictides has been discussed 
previously.^ For the different behavior of the Cg and t2g 
orbitals it is important that the Cg states in LiFeAs lie en- 
ergetically lower than the t2g states. In the atomic limit, 
the energy gain from Hund's rule exceeds the crystal field 
splitting already for rather small J and the ground state 
is a high spin state with the configuration egt2g where the 
t2g orbitals are occupied by three electrons of the same 
spin. In the atomic limit this prevents mixing of the or- 
bitals due to the Pauli principle; in the crystal, it still 
impedes inter-orbital fluctuations within the t2g mani- 
foldj ^^ l ^^ ^This effect contributes to the high sensitivity of 
the t2g effective masses with respect to J. 
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FIG. 8. (Color online) Sensitivity of the Fermi surface at 
kz = with respect to changes in the interaction parameters. 

Accordingly, the Fermi surface is rather stable against 
variation of U but strongly depends on J, shown for 
kz = in Fig. [S] and A:;^ = tt in Fig. [9] Following the 
trend discussed above, larger values of J promote a more 
pronounced shrinking (increase) of the middle (outer) 
hole pocket. Values as large as 0.9 eV for the Hund's 
coupling render the system rather incoherent though, 
with significant scattering rates -ImE(zO+) around 14 
meV on the t2g orbitals. This causes the broadening 
of the respective Fermi surface in Figs. [S] and |9] Future 
calculations with the full Hund's rule coupling are re- 
quired to check whether the occurrence of the coherence- 
incoherence crossover already at J < 0.9 eV is a physical 
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U=3.5 eV, J=0.8 eV U=4.5 eV, J=0.8 eV 
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FIG. 9. (Color online) Sensitivity of the Fermi surface at 
kz = TV with respect to changes in the interaction parameters. 

effect or an artifact of the breaking of the rotational in- 



variance by density-density interactions. 

IV. CONCLUSIONS 

We showed that within the considered range of in- 
teraction parameters LiFeAs behaves as a Fermi liquid 
where correlation effects are very sensitive to the value 
of the Hund's rule coupling in particular for the t2g or- 
bitals. The strong mass enhancements measured in both 
ARPES and dHvA experiments suggest sizable correla- 
tions of the size considered in this work. While electron- 
phonon effects have been reported^^ to be significant and 
contribute to the slightly higher mass enhancements mea- 
sured experimentally, they cannot account alone for the 
large values observed. As for the Fermi surface, the cor- 
relations mainly affect the hole pockets that significantly 
change in size. We propose this as the source of the seem- 
ing discrepancy of the ARPES and dHvA experiments: 
whereas dHvA presumably observes only electron orbits 
with sizes close to their LDA values, ARPES finds the re- 
duced size of the middle hole pocket as most prominent 
feature. In this way, the two experiments can be rec- 
onciled. The selective size reduction of the middle hole 
pocket also renders nesting less efficient. 
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